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Abstract 

We consider a quantum system constituted by N identical particles interacting by means of a mean-field 
Hamiltonian. It is well known that, in the limit N —^ oo, the one-particle state obeys to the Hartree equation. 
O , Moreover, propagation of chaos holds. In this paper, we take care of the h dependence by considering the 
semiclassical expansion of the A^-particle system. We prove that each term of the expansion agrees, in the limit 
OO ' N — > cxD, with the corresponding one associated with the Hartree equation. We work in the classical phase space 
by using the Wigner formalism, which seems to be the most appropriate for the present problem. 
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^ ■ 1. Introduction 

The Hartree equation is the following nonlinear one-particle Schrodinger equation: 

^n9i^ = -yAV^+(0*|^l')V^, (1.1) 

where the mass of the particle is chosen equal to one, : ]R3 ^ M is the two-body potential 
and 

{<P*\^\'){x)= [ dy<P{x-y)my)\' (1.2) 
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is the effective self-consistent interaction. 

Equation fll.ip arises as the reduced description of a system of identical particles interacting 
by means of the mean-field potential: 

1 ^ 

U,,{X^) = - J2 <f>i^i-^^)^ (1-3) 

l<«<jr<Af 

(where = {xi, . . . , xn}, Xj G M^, j = 1, . . . , A^) in the limit N ^ oo. 

In fact, consider the A^-particle wave function = \l/Ar(XAr;t) solution of the Schrodinger 
equation: 

1=1 

with a completely factorized initial state given by: 

N 

^n{Xn; 0) = ^n,o{Xn) = llM^j)- (1-5) 

Then, it is well known that the j-particle reduced density matrices, defined as: 

pf (X„ Yf, t)= [ dX^_,^jv iX„ X^-r, t) {Y„ Xn-f, t) , (1.6) 

7lR3(JV-j) 

converges, in the limit N ^ oo and for any fixed j = 1, . . . , A^, to the factorized state: 

i 

k=l 

where il){x; t) solves the one-particle Hartree equation (11. ip with initial datum ipQ. This feature 
is usually called "propagation of chaos". 

The previous result was originally obtained for sufficiently smooth potentials (see [1], [2], 
[3]); then it has been generalized to include Coulomb interactions (see [5], [6]). Furthermore, 
some results concerning the speed of convergence of the mean-field evolution to the Hartree 
dynamics (for all fixed times), have been proven more recently (see [7j, |8j). 

The limit N —>■ oo for a classical system interacting by means of the same mean-field 
interaction (11. 3p . can be considered as well (see 0, [ID], [H], [12], [13] for the case of smooth 
potential, and [14J for more singular interaction). In fact, considering as initial state of the 
system a probability measure F/v o = -F/v,o {X^, Vn) in the A^-particle phase space which is 
completely factorized, namely: 

N 

Xj, Vjj , l-L-oj 



where /o is a given one-particle probability density, it is known that its evolution at time t > 0, 
denoted by F^- {X^, V^; t), is obtained by solving the Liouville equation: 

(at + V)vVx^)Fjv = Vx^f/iv-Vy^Fjv. (1.9) 
Then, the j-particle distribution at time t > 0, defined as 

/f (X„ V,; t)= I dXjv_,d\/^_,F^ (X„ V,, Xn_„ Vn-,; t) , (1.10) 

converges, as oo, to the product state: 

j 

f,{X„V,-t) = l[f{xk,Vk;t), (1.11) 

fc=i 

where f{x,v; t) is the solution of the Vlasov equation: 

{dt + v-V,)f = {V4>*f)-VJ, (1.12) 

(the convolution above is with respect to both the variables x and v) with initial datum /q. 
Equation f ll.l2p is the classical analogous of the Hartree equation (11. ip . 

Although the mean-field limit ^ oo is well understood for both classical and quantum 
systems, there is a question which seems to be still open, namely, does that limit hold for 
quantum systems uniformly in h, at least for systems having a reasonable classical analogue? 
The proofs which are available up to now exhibit an error vanishing when N —>■ oo but diverging 
as ^ — > 0, although in [16], [17], [H], [19] some efforts in the direction of a better control of the 
error term have been done. 

If one wants to deal with the classical and quantum case simultaneously, it is natural to 
work in the classical phase space by using the Wigner formalism. 

The one-particle Wigner function associated with the wave function iIj{x; t) is given by: 

f (x, v; t) = (27r)-3 £^ dy e^^-> (^x + ^;t^^p(^x- ^; , (1.13) 

and, similarly, the A^-particle Wigner function associated with the wave function N{XN;t) is 
defined as: 

(X^, VW; t) = (27r)^3^ 1^^^ dYj, e^^-^-^^ + " ^) • (l-^^) 

Then, by using that ip{x; t) and \E'Ar(X7v; t) solve equations (11.10 and (II. 4p respectively, we find 
the equations: 

{dt + vV,)f = T^f (1.15) 

and 

{dt + Vn ■ Vx^) = T^W^, (1.16) 
where and are suitable pseudodifferential operators. 



The initial data for equations f ll.lSp and (11.161) are 

(x, v) = (27r)-=^ £ dy (^x + ^ j - y) , (1-17) 



and 



N 

respectively. 

One can easily rephrase the result of [I] by showing that the j-particle Wigner function 
W^^ (X,, Vf, t)= f dXN-,dVN-jW^ {Xj,Xn.,, V,, VN-f, t) (1.19) 

converges, in a suitable sense, to 

j 

ff (X„ Vj- t) = \[f (xfc, Vk] t) for any t > 0. (1.20) 

k=l 

However, the convergence error is diverging when — > (for example, for sufficiently small 

■ c hp) 

times t < to it is of the form -^efi). The reason is that the operators T and appearing 
in f ll.lSp and (11.161) are bounded as operators acting on the space in which we can prove the 
convergence of (ll.lQp . but their norm diverge as | when /?, — > 0. On the other hand, the classical 
counterpart of this problem has been solved, so that it seems natural to look for an asymptotic 
expansion for the j-particle distributions j, namely: 

it) = <i (t) + (t) + h'W^^] (t) + . . . , (1.21) 

and for an analogous expansion for the j-fold product of solutions of the equation (ll.lSp . 
namely: 

f- it) = ifT' (^) = (t) + ^ff (^) + ^'ff (^) + • • • (1-22) 
The zero order term in (ll.2ip corresponds properly to what we previously denoted by /j^, while 

the function /j"'' appearing in (ll.22p is exactly the j-fold product of what we called / (see (ll.lOp 
and (I1.12P ). Therefore, at zero order in h we obtained the classical quantities, as expected, and 
we know that the convergence of W^"*- to f^^^ is well established. Then, it looks natural to try 
to show the convergence 

W^^]it) ^ ff\t) , as ^ oo, for any A; > 0. (1.23) 

This is the goal of the present paper. 

A complete proof of the uniformity in h of the limit N oo would require a control of the 
remainder of the expansion (11.210 . but we are not able to do this. However the error term of 
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order k in the expansion f ll.22p can be proven to be 0{h''^^) by adapting the proof in [15] for 
the hnear case to the present context, under suitable smoothness assumptions. 

The plan of the paper is the following. In the next two sections we present a semiclassical 
expansion of the Hartree equation and of the iV-particle system. After a brief discussion of 
the hierarchical structures and their inadequacy as a technical tool for the present problem, 
we introduce the classical mean- field limit which is the basis of our analysis. After that, we 
explain the strategy of the proof of the convergence (I1.23P . The last two sections are devoted to 
the precise statement of our result, its proof and supplementary comments. Three Appendices 
contain technicalities. 



2. The Hartree dynamics 
The Wigner-Liouville equation associated with the Hartree equation (11.11) reads as: 

(9, + .-v..)r = T;.A 

where, from now on, we set e = h. Furthermore {x, f ) G M'^ x and /^(x, v; t) is defined as in 
( I1.13P and it is a real (but non necessarily positive) function on the classical phase space. We 
remark that in the previous section we used the notation for the right hand side of (12. ip 
but now we use T^^f^ to stress the nonlinearity of equation (12. ip . 
Here acts as follows: 

(x, v) = {2ti)'H f'^ dA / dk^{k)pg{k)e ''■^{k ■ V,)r (x, V + eXk), (2.2) 

J-l/2 Jr3 

where, as usual, = 0(x) denotes the two-body interaction potential and 



Pg{x) = dv g{x,v). (2.3) 

Finally, we denoted by Pg and the Fourier transforms of Pg and (f) respectively, namely: 

Pg{k) = [ dx e-' ^""pgix) and 0(A;) = [ dx e"^ '^■^0(x). (2.4) 

JR3 J^3 

Following p!3], for a fixed g, can be expanded as 

T; = Tf + £T^(i) + e^T^ + ... (2.5) 



where 



= Cr,{2n)-H [ dkmPgWe' '-^k ■ (2-6) 

(2.7) 



" 2"(n + l)!' 
for n even and 

T(") = 0, (2 



for n odd. The operator Tg^\ for n even, can be also written as 

) = (-i)"/2c„ (D^V * g) ■ d:^\ (2.9) 

where, as in (11.121) . * denotes the convolution with respect to both x and v and we used the 
notation: 



D'^vD''C= V — — (2 10) 



ni,n2,n3£N: 
T,jnj=n 

with X = (x\x^x^) e and v = {v^,v^,v^) G 

for the one-particle functions u and (. 

We want to determine a semiclassical expansion of /^(x,f;t), solution of equation (12. ip . 
namely: 

fit) = f<'\t) + £/«(t) + e'f'\t) + ... (2.11) 

First of all, we observe that interesting physical states, such as coherent states (see e.g. the 
forthcoming Section 6), have a similar expansion even at time t = 0, namely: 

/o^ = /r+^/o^'^+^Vo^'^ + --- (2.12) 
Inserting ( 12. lip in ( 12. 5p and setting: 

Tt^ = Tit (2.13) 

we readily arrive to the following sequence of problems for the coefficients f^''\t) of the expan- 
sion (121T1) : 

(9, + ■ V.) /(°) = t( 
f^'\x,v;t)l^^ = f!,'\x,v), 

and 

(5t + -v.) = L(/W)/('^) + e('=), 

f^''\x,v;t)l^^ = fl,'\x,v), 
for k > 1, where 

L{h)f = Tl^^f + Tfh = (V,0 * h) ■ VJ + (V,0 * /) ■ V,h, (2.16) 

and 

_ ^ T,(^)/«. (2.17) 

r,s,l: 
s+r+l=k 
l<k,r<k 

Note that equation fl2.14p is nothing else than the classical Vlasov equation (see (11.120 ) associ- 
ated with the interaction (p (which will assume to be smooth). It is well known that the Vlasov 
equation can be solved by means of characteristics and fixed point. Moreover, the problems 
(I2.15P are linear and can be solved by a recursive argument (see Section 6 below). Indeed, the 
source terms Q^''^ involve only those coefficients /'•"^ with n < k, so that they are known by 



(2.14) 



(2.15) 



the previous steps. We shall give a sense to the solutions f^''\t), for /c > 1, once we will have 
established precise assumptions on and /g'^'*. 

3. The A^-particle dynamics 

Consider a quantum system constituted by N identical particles. Its time evolution in the 
classical phase space is given by the Wigner-Liouville equation: 

{dt + Vn ■ Vx^) Wl, = TI^WI,, (3.1) 

where W^{Xn, Vn', t) is the Wigner function describing the state of the system, 

Xn = {xi,...,xn) ew'"', VN = {vi,...,VN)em.^'', 

and the pair Zjs; := (X^?, V^v) denotes the generic point in the classical A^-particle phase space. 
Moreover, 



(3.2) 

where K]\r = {ki, . . . , k^) G and 11^ is the (mean-field) interaction potential given by: 



N 



l<l<j<N 

We choose, as initial datum, the factorized state: 

N 

W^^^,{XM,VM) = l[f^{x„v,), (3.4) 

where /g is the initial datum of the equation fl2.ip . Following [T3], we expand 

= + eTl}'> + e'T^ + ... (3.5) 

where, for n even we have 

4") = z(27r)-3^C„ / dir^f/(ir^)e^^--^- {K^ ■ Vy^T^' , (3.6) 

Cn being constants depending on n, and, for n odd, we find 

T^"^ = 0. (3.7) 

Looking for a semiclassical expansion 

W'^it) = wPit) + eW^^\t) + e'wg\t) + ..., (3.8) 
we first expand the initial datum 

W'n,o = <o + + ^'<o + • • • • (3.9) 
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r{k) 



The coefficients W^jq are determined by (13 .4^ and (I2.12p as 



N 



W^Z = \{fo\^3^v,l (3.10) 



AT 

<i= E n/o'^^(^^'^^)- (3-11) 



0<Sj<fc 

Note that w'^\ is factorized only for A; = 0. 

By (13. 8p and (13. Sp . we arrive to the sequence of problems: 



(3.12) 



t=o 



and 



(9, + ■ Vx.) w^J:^ = + e^'^) 



N ' 



wj,'\XM, Vn; t) = W'^'>,{Xn, Vn) 



(3.13) 



i=0 



for k > 1, where 



^(fc) _ ^ ^(fc-OT;f^(Z) 

0<«<fc 

.(0) 



©iJ^ = E (3.14) 



Note that = VxmUn ■ Vy^ = ;^ X^kkjxa^ ^2;0(a;i — a;^) ■ Vi,. is the classical Liouville 



operator, while the source terms O^'*, at each order k, are known by the previous steps. Note 
also that, under reasonable assumptions on the interaction potential 0, equation (13.120 can be 
solved by considering the Hamiltonian flow $*(X7v, V/v), solution of the problem 



N 



~^^x(p{.Xi-Xj] 



(3.15) 



Indeed 

wP{Xn, Vn- t) = SN{t)wP^^{XN, Vn) = W^^l ($-* (X;v, Vn)) , (3.16) 

where, from now on, we denote by Sn the flow generated by the Liouville operator while 
equations (13.130 can be solved by recurrence thanks to the Duhamel formula: 

ft 



W\;\t) = 5^(t)W^jV,o + / dti SNit - h)Q';>ih). (3.17) 







We conclude this section by expressing the operators T^"* [n even) in terms of the variables 
Xn, Vn- From fl3.6p . we find that: 

4") = + i?Sv"\ (3.18) 

where 

/I \n/2 

T^N^=Cn^-^ E ^rV(a:.-a:,)-/^.r\ (3.19) 

l<i<j<iV 

where c„ is the same of (12. 7p . and 



b: 



( ) 1 - Qn+l Qn+l 



l<l<j<N fci,fc2eN3 ^xi uxj 

\ki\ + \k2\=n+l 



where, for i = 1,2, ki = {ki^i, ki^2, ^3), \ki\ = Ki + + ki,3 , and 

gn+l g\ki\ g\k2\ 



(3.21; 



I "^j L L t J J J 

with Xj = [x], x^, x^j) , xi = [x], x], xf) , 

while Cki^k2 are suitable coefficients. The same holds for the derivatives with respect to the 
velocities. 

We observe that, by the expression (13.201) . we mean that the derivative of order \ki\ is 
distributed over the three components of xi in the same way in which it is distributed over the 
three components of vi, and the same holds for the derivative of order \k2\- 

4. Hierarchies 

One way to investigate the behavior of the A^-particle system in the limit — > 00, is to 
consider the hierarchy associated with equation (I3.ip . More precisely, introducing the j-particle 
functions: 

Ty^^.(X„ Vf, t)= ! dXN-,dVj,.,W%{X„ V„ V^_,; t), (4.1) 

a straightforward computation yields the following sequence of equations: 

[dt + V, ■ Vx,) W^^ = T^WI,^^ + (^^) q+iW^^,,+i, J = 1, 2, . . . , AT, 



(4.2) 



with ly^^AT = Wf^ and C^+i = 0, 



which, as usual, is called "hierarchy" because each equation is linked to the subsequent one. 
The operator TJ (for a fixed j) describes the interaction of the first j particles (we recall that 
we are dealing with identical particles, thus, in order to refer to any group of j particles we can 
say "the first j particles" because we can rearrange them as we want) among themselves, while 
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the operator describes the interaction of the first j particles with the remaining N — j. 
The exphcit form of such operators is: 



y [ dX [ dk4>{k)e'''<^^-^^\k-V.,)W'^JX,,Vi.i,vi + Xek,Vj.i), 
-TT^^J -1/2 Jr3 

(4.3) 



and 



i 1-1/2 



ik-{xi-Xj+i) 



= z(27r)-=^^V / dA / dk4){k) [ dXj+idVj 

{k ■ V^,i)W^j^-i^{Xj,Xj+i,Vi-i,vi + Xek,Vj-i,Vj+i). 

(4.4) 

Note now that Tj = O (^^j; thus we expect its action to be neghgible in the hmit. On the 
other hand, if we consider the sequence {fj(t)}j, where fj(t) = fj{Xj, Vj; t) is given by: 

j 

f^{X„V,;t) = l[nxk,Vk;t) (4.5) 

k=l 

and /^(t) is the solution of the Wigner-Liouville equation associated with the Hartree dynamics 
(see (12.11) ), we easily deduce the following hierarchy: 

{dt + VrVx,)f] = q^J]^,. (4.6) 

Therefore, we expect that 

asiV^oo, (4.7) 

for any t > 0, provided that the same convergence holds at time t = 0. As we have already 
recalled in Section 1, such a result can indeed be proven under various assumptions on the 
interaction potential 0, both in the reduced density matrix formalism and in the Wigner one. 
This constitutes a validation of the Hartree equation in the mean-field limit. Nevertheless, as 
we have already remarked, a common features of these results is that the limit is singularly 
behaving when e —>■ 0. In fact, the operators involved in the above hierarchies are bounded in 
the norm appropriate to study the convergence, both in the reduced density matrix formalism 
and in the Wigner one, but their norm is diverging when e goes to zero. This suggests to consider 
the semiclassical equations described in Sections 1 and 2. In this way, considering equations 
at each order in e and analyzing the hierarchies associated with each of those equation, we 
have to deal with operators which are clearly independent of e {T^\ for the iV-particle case, 
and T^"''' for the expansion associated with the Hartree dynamics), and we have to investigate 
only the limit N oo without any dependence on e. The price we have to pay is that now 
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those operators are unbounded, as it comes out for the classical mean-field limit we are going 
to discuss in the next section. Therefore, if we want to prove that the coefficient of order e'' of 
the expansion of VV^j, namely: 

W^^'UX,, V,- 1)= ! dX;v-,dVW-,W^j?^(X„ X;v-„ V,, VM-y, t), 

jR3{]V-j)xR3{]V-j) 

(4.8) 

converges to the corresponding object for the Hartree fiow, which is: 

fj'\x,,Vr,t) = J2 tlf^'^\xr,vr;t), (4.9) 

si...Sj: r=l 

0<Sr<k 

(where the one-particle functions f^^''\xr,Vr;t) solve equation (12.141) . if Sr = 0, and fl2.15p . for 
Sr > 0), the use of the hierarchy solved by Vr^](t) does not seem a good idea. In fact, even 
at level zero, when we have to deal with the classical mean- field limit, the hierarchy is very 
difficult to handle with. Such a hierarchy is given by: 

(S, + V, . V.J <. = 7;'°'<. + ^cfA<.«, (4.10) 
(i;'°'<;)(.Y„l/,) = l Yl V,,,#.(i, - a:,) • V„M/<,°>(A',, VJ, (4.11) 

l</<r<j 

(cfil<J,.)(A,.l{,) = 

Clearly these operators are unbounded (unless to make them act on analytical functions) be- 
cause they involve derivatives with respect to the velocity variables. For this reason, to deal 
with the hierarchy is quite difficult and the obstacle which occurs in facing the higher order 
terms is precisely the same. However, in the classical case we can treat the convergence in a 
more natural way, avoiding to use the hierarchy. The idea is to control the j-particle marginals 
W^^j in terms of the expectation of the j-fold product of empirical measures with respect to the 
initial iV-particle probability distribution (see Section 5 below). In the present paper we follow 
a similar strategy in dealing with the convergence of the higher order terms of the expansion. 
More precisely, we will express W^'^j in terms of the expectation, with respect to the initial 
iV-particle zero order coefficient (which is known to be a probability distribution), of suitable 
operators acting on empirical measures. The control of these objects will be obtained thanks 



with 



and 



12 



to some estimates of the derivatives of the classical flow with respect to the initial data (see 
Proposition 5.1). 

5. The classical mean-field limit 

The semiclassical expansion of the A^-particle system leads us to consider the sequence of 
problems (13. 121) - (13.131) . The zero order equation (I3.12p is purely classical and well understood. 
At this regard, let us remind some basic facts concerning the case of smooth potentials, which 
will be crucial in what follows. 

The Vlasov equation (11.121) makes sense even for a generic probability measure v because 
V0 * 1/ G C°°(M^) (thanks to the smoothness of 0) and, by using the characteristic flow 
associated with the equation and a fixed point argument, it is possible to prove the existence 
and uniqueness of the solution. Furthermore, introducing the Wasserstein distance W (e.g. 
[T3] ) based on a bounded metric in to avoid unnecessary boundedness assumptions on the 
moments of the measures we deal with, it is possible to prove the following continuity property: 

W{vlvj)<e'''W{vlvl) (5.1) 

where z/q and Vq are two probability measures and v} and v"^ are the weak solutions of the 
Vlasov equation with initial data given by z/q and z/q respectively (see again [I3j). Moreover, 
for a configuration = {zi, . . . , zn}, where Zj = {xj,Vj) G M^, consider the Hamiltonian flow 

<l>*(X;v, Vn) = Z^{t) = Z^it- Zn), (5.2) 

with initial datum Z^ (see (13.151) ). and construct the empirical measure ^Nif) as follows: 

1 ^ 

/i7v(t) := ^iN{z\ZN{t)) = j;^^^{z- Zj{t)) . (5.3) 

j=i 

The basic remark is that fiN{t) is a weak solution of the Vlasov equation, so that, by (15. ip . 
we have: 

Mt)^f^''\t), asiV^oo, (5.4) 

provided that 

/i^ - /^°H0) = /o^°\ asiV^oo, (5.5) 

where 

1 ^ 

fJ^N ■= fJ'N{z\ZN) = {Z - Zj) (5.6) 

i=i 

is the empirical distribution at time t = 0. Moreover, Jq^^ is a (possibly smooth) probability 
distribution and f^^\t) is the solution of the Vlasov equation with initial datum /q^'*. Clearly, 
the convergences (15. 4p and (15.50 hold with respect to the metric induced by W on the space of 
probability measures on and this is equivalent to the weak topology of probability measures. 
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Next, let us consider the (factorized) A^-particle probability distribution W^^qIZn] 



Ylk=i fo i^k, Vk) and let Wj^ {t) be its time evolution according to the Liouville equation fl3.16p . 
We want to investigate the behavior of the j-particle marginals W^j{t). Denoting by E^v the 
expectation with respect to wj^\{ZN), after straightforward computations, we obtain: 

[l^N iz[\Z^it)) . . ./i^. (4|Z^(t))] = NiN-l)..^iN~j + l) ^(0)^^..^) + 

(5.7) 

where = {z[... z'^) and w''^]{Z'r t) = wj^\{Z'^{t)) (see 1^). 

Consider now a typical sequence Z^ with respect to /q°\ namely such that (15. 5p holds. By 
the strong law of large numbers this happens a.e. with respect to ( /o ) and by (15.41) and 
(15.71) we have: 

hm E^ [fx^{z[\Z^it))...ix^ (4|Z^(t))] = hm W^^]{Z'.,t) = {f^'^f (Z'-^t), 

(5.8) 

in the weak topology of probability measures. Thus propagation of chaos is proven, and, this 
is the remarkable fact, it has been done without using the hierarchy. 

For fixed e > 0, the quantum hierarchy is, in a certain sense, easier. In fact, in that 
situation the operators involved are bounded (as operators acting on the spaces appropriate for 
that context) and it is possible to realize the limit by using the hierarchy. On the contrary, in 
the quantum context we cannot use any characteristic flow and there is not any object analogous 
to the empirical measure. Nevertheless, if we consider the semiclassical expansion of the time 
evolved Wigner function, the higher order terms can be viewed as quantum corrections to the 
classical dynamics. 

We now explain heuristically our approach fully exploited in Section 7. 

The first correction to the Vlasov equation in the Hartree dynamics satisfies (see (12.111) 
and fl2:T2|) ): 

(9, + .-V.)/« = L(/(°))/«, 

(looking at the expression (I2.17P for the source terms G^^\ we verify that G^^-* =0). As we shall 
see in detail in the following section, our choice for the initial one-particle datum is a mixture 
of coherent states and each coefficient of the expansion it is given by suitable derivatives of the 
zero order distribution. In particular, the explicit form for /q^^ is: 

f!,'\x,v) = Dyl,'\x,v), (5.10) 
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where Dq is a suitable second order derivation operator (see formula (16.81) below in the case 
k = 2) involving derivatives with respect to the initial variables zi, . . . , z^- 

As regard to the A^-particle dynamics, looking at (13. lip in the case k = 1, we know that 
the initial datum for the coefficient of order one in e is: 



N N 



<o(^^) = E /o'^(^^) n /o°^(^') = ^'<o(^^)' (5-11) 

where 

N 

^' = E^G,,' (5.12) 
i=i 

and Dqj is the operator Dq relative to the variable Zj G MP . Let us now define V^fiN^t) as 
the distribution acting on a test function u in the following way: 

(n,PV^(t)) = V' U$^^(^Kt)) = ]^ E ^lA'iit))- (5-13) 

\ 1=1 J 1,3=1 

We remind that the operator D'q ^ involves derivatives with respect to the initial variables 

^1, . . . ,2;7V, thus, if at time t = we have f^^ when ^ oo (in the weak sense of 

probability measures), it follows that: 

1 ^ 1 ^ 1 ^ 

1=1 i,j=i j=i 

(5.14) 

as ^ oo. Moreover, by (15. lip and (I5.14p . we can conclude that: 

(5.15) 



By equation (13.130 for k = 1, we have: 



w'}^\ZM-^)\t=, = W^^I{Zm), (5.16) 
namely, the classical Liouville equation. Therefore: 

W^^\ZN]t) = SN{t)W^^\{Zx). (5.17) 
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Finally, by virtue of f l5.17p and f lS.lip . we obtain 



' JR3]VxM3JV 

R3^xK3JV 



R3JVx 



(Zn) {u,fiNit)) 



R3iVxM3JV 

= {u,EN[V^fiNit)]). (5.18) 

Therefore, the behavior of wl^\{t) is determined by that of T>'^ijLj^{t) for any initial configuration 
Zjv which is typical with respect to fj^\ Finally, since /iAr(t) solves: 

(5.19) 



applying "D^, we get: 



'^'lM\LI=ljiaAHbllL^ UMiLI -t- n.M. 

(5.20) 



{dt + v- V.) V^fLNit) = L (/i^(t)) PViv(t) + Rn 



where Rn is a term involving objects of the form {DcjfiNit)) {DcjUNif)) which, as we 
shall see later, are of order 1/N when tested versus smooth functions. The equation fl5.20p is 
similar to (15. 9p . except for the presence of the term R^ and for the fact that we have L {^iy(t)) 
instead of L [f^^^). Therefore, the proof of the convergence of W^^\{t) to /(i)(t) reduces to 
that of a stability property for the solution of (15. 9p with respect to suitable weak topologies. 
Propositions 6.1 and 6.2 below will provide us such property. 

The general case /c > 1 is only technically more complicated because of the presence of 
source terms, but the main ideas are those presented here. 

We conclude by establishing a Proposition controlling the size of the derivatives of the 
Hamilton flow (I3.15P with respect to the initial data. 

From now on we shall denote by C a positive constant, independent of N, possibly chang- 
ing from line to line. 

Proposition 5.1 

Let Zi{t) = {xi{t),Vi{t)) , i = 1,...,N be the solution of equations Ii3.15\) with initial datum 
Zi = (xj, Vi) , i = 1, . . . , N . Let zf W [3 = 1, . . . ,6 be the f3-th component of Zi G M^. // the pair 
interaction potential is C°°(R'^) and the derivatives of any order of (p o^^e uniformly bounded, 
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then, for each k G N; 



< 



c 



(i) 5 



(5.21) 



where I := is any sequence of possibly repeated indices and d^^ is the number of 

different indices in I which are also different from i. 

The physical significance of (15.211) is obvious. In the mean-field context, the quantity Zi{t) 



depends weakly on Zj ii j ^ i for each t > 0. Actually 



O (^) while 



Oil) and 



J » 

these two estimates give rise to (15.211) in the case k = 1. Estimate (I5.2ip says that for each 
derivative of any order with respect to some Zj of Zi{t) , we gain a factor 1/A^. We have also 
the following corollary whose straightforward proof will be omitted. 

Proposition 5.2 

Let U = U{ZN{t)) be a function of the time evolved configuration Z^it) of the form: 

N 



[/(Z^(t)) = lX^n(z.(t)) 

i=l 



where u E C2 



X 



) . Then, if the pair interaction potential satisfies the assumptions of 



Proposition 5.1, the following estimate holds: 

d^u{ZM{t)) 



UZj^ . . . UZj^ 



< 



C 



(5.22) 



where dk is the number of different indices in the sequence I = (ji, . . . ,jk)- 
The proof of Proposition 5.1 will be given in Appendix A. 



6. Results and technical preliminaries 

We choose, as initial condition for the one-particle Wigner function, a mixture of coherent 
states. The Wigner function associated with a pure coherent state centered at the point (xq, Vq) 
is given by: 

w{x,v\xo,vo) = — r^e e e = . (6.1) 

Let now g = g{x,v) be a smooth probability density on the one-particle phase space 
independent of e (see Hypotheses H below) . Then we define: 

fo{x,v) = / dxodvo w{x,v\xo,vo)g{xo,vo). (6.2) 
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Using the standard notation z = {x,v) and zq = (xo,fo), (I6.2p is equivalent to: 

1 



[ dCe-^'g{z-V^O. (6.3) 



(TT 

Expanding 

- v^C) = 9i^) - (C ■ V.) giz)V~e + (C ■ V.fgiz)^-^ + ... 

(2n-l)! ' (2n) 



(6.4) 

and performing the gaussian integrations (which cancels the terms with the odd powers of \^), 
we readily arrive to the following expansion for the Wigner function /q: 

/o = + efl,'^ + ■■■ + eVt^ + ..., (6.5) 

where 

/o^°^ = 9, (6.6) 

= Dl^fS"^ for n > 1, (6.7) 

and Dq {G stands for "Gaussian"), for each /c > 0, is the following derivation operator with 
respect to the variable z = {x,v): 

(9z"i . . . 

aj=l,...fi 

where 

k 

CGiai...a,) = ^ [ dCe-^'TTO. (6.9) 

Therefore, Ccio^i ■ ■ ■ ctfe) is equal to zero for each sequence ai . . . in which at least one index 
appears an odd number of times. 

Hypotheses H: 

In the present paper we assume that the probability density g = f^^^ G S{M.^ x R^), thus 
(16.71) make sense for any n > 1. As regard to the pair interaction potential 0, we assume 
that (j) G C°°(R^), that any derivative of is uniformly bounded (in order to be able to apply 
Proposition 5.1) and that is spherically symmetric. 



D^^= V CG(ai...afc)- , (6.^ 
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Remark 6.1: 

In this paper we consider a completely factorized A^-particle initial state. Furthermore the 
one-particle state is a mixture and this automatically excludes the Bose statistics. 

Remark 6.2: 

We made the choice to expand fully the initial state /o according to equation (16.51) . An- 
other possibility is to assume the {e dependent) state /o (which is a probability measure in the 
present case) as initial condition for the Vlasov problem and, consequently, /o = for the 
problems (I2.15p . Now the coefficients f^'^^t) are e dependent but this does not change deeply 
our analysis because /o is smooth, uniformly in e. 

Under hypotheses H, we can give a sense to the linear problem (12.151) for any /c > 1, by 
virtue of the following proposition, whose (straightforward) proof will be given in Appendix B. 

Proposition 6.1 

Consider the following initial value problem: 



with 7o G L1(M3 X R^), h = h{x,v;t) is such that \V^h\ G C° {L\M.^ x M3)^k+), 6 = e{x,v;t) 
is such that 6 G C° (^^(R^ x M3)^k+). 

Then, there exists a unique solution'-^ = 'j{x,v;t) of Ii6.10\) . such that'-/ G {L^{M.^ x ]R^),R+), 
given by an explicit series expansion. 

Furthermore, denoting by the flow generated by L{h), we have that T,h(t, 0)70 G C'^ {M.^ x M.^) 
provided that V^h G (M^ x M^) and 70 G (M^ x R^). 

The main goal of the present paper is to compare the j'-particle semiclassical expansion 
associated with the A^-particle flow, namely W^j{t), k = 0,1,2,..., with the corresponding 
coefficients f\'^\t) of the expansion: 



Theorem 6.1 

Under the Hypotheses H, for all t > 0, for any integers k and j , the following limit holds in 
S'(R^^ X M^i). 




(6.10) 



= f^'\t) + effit) + ■■■ + e'ffit) + . . . , 



(6.11) 



where /j (t) is given by (14.91) . The main result is the following. 



(6.12) 
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as iV — >• oo. 
Remark 6.3: 

As we shall see in the sequel, the convergence (16.121) is slightly stronger than the con- 
vergence in iS'(M^-' X M^-'). Indeed, the sequence converges also when it is tested on 
functions in C^{M.^^ x M.^^), namely, the space of functions which are uniformly bounded and 
infinitely differentiable. Such kind of convergence, which is natural in the present context, will 
be called C^-weak convergence. 

Proposition 6.2 

Let 7Ar(a;, be a sequence in S'(M.^ x R'^) (for each t) satisfying: 

{dt + v ■ Vx) In = L{hN)jN + ©tv, , , 

/ \i / \ (6.13) 

77v(a;,f;t)li=o = 1n,o{x,v), 

where 7Ar,0) ©iv sjre sequences in iS'(]R^ x R^). We assume that: 

i) hp^{x,v;t) is a sequence of probability measures converging, as N ^ oo, to a measure 
h{t)dxdv with a density h{t) e (R^ x R^) and such that \V^h\ e C° (L^(R^ x R3),R+). 

ii) for all Ui,U2 in C^(R^ x R^) , there exists a constant C = C{ui,U2) > 0, not depending on 
N, such that: 

* (u277v)||ioo(K3xR3) < C < +CX) for any t. (6.14) 

iii) jn,o 70) as N ^ oo, C^-weakly , 70 = 7o(x, f) is a function in L^{R^ x R^). 

iv) Qn -^Q, as N ^ 00, C^-weakly , © = Q{x,v;t) is a function m C^{L^{W^ x R3),R+). 
Then: 



In ^ 1, as N 00 -weakly, 

where 7 is the unique solution of the problem Ii6.10\) in C° {L^ (R'^ x R'^) 
For the proof, see Appendix B. 

7. Convergence 
This section is devoted to the proof of Theorem 6.1. 

By flaTTD and (EJl), for A; > we have: 

n>0 r=0 n...r„: 
r,>0 



(6.15) 



Jo Jo 



rj=k—r 

SN{t - t,)Tj^'^SN{ti - t2) . ..T'^;-^SN{tn)W^;;\{ZN). (7.1) 
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It is useful to remind that the only non-vanishing terms in (17. ip are those for which all 
ri, . . . , r„ are even. 

According to (13. lip and (16. 7p . 

<o(^iv)= E fl{D'c:4"\'^))^ (7-2) 

SI...SN 7 = 1 

0<Sj<r 

where Dqj is defined in (16. 8p and the extra symbol j means that this operator acts on the 
variable Zj. Defining the operator as: 



N 



Si...sjv: = 1 
0<Sj<r 



we have: 



iy;7o(Z^) = P^W;^^(Z^) V r > 0. (7.4) 

In order to investigate the behavior of the j-particle functions W^''-{Zj]t) when N — > oo, 
we consider the following object, for a given configuration = (2;^ . . . 2;^): 

a;}f^^.(Z;.;t)= / dZ^ W^^\ZN■,t)^^M{z[\Zr,) . . . ^lN{z'^\ZN). (7.5) 

JR6JV 

In the end of the section, we will show that (17. 5p is asymptotically equivalent to wl^j{Zj] t). 
From ([71]), ([72D and ([73D, it follows that: 

ijP^iZ'^; t) = y^y^ V dti r dh... r ' dtn [ dZMiiNAZ'^\ZN) 



n>0 r=0 

r,>0 
E rj=k—r 

SM{t - h)Tj^'^SN{h - h) . . .4^"^^^(t„)p2wi°) (Z^), (7.6) 

where 

^iN,j{Z'AZN) = hn{z[\Zn) . ..hn{zj\Zn). (7.7) 
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(r ■) 

Integrating by parts, reminding that each vj is even and that each involves derivatives of 
order rj + 1, we have: 



n>0 — " " ■ - ^ord 



E 



|r„|=fe-r 



where r„ is the sequence of positive integers ri, . . . , r„, |r„| = X]j=i '"j ZAr(t) is the Hamil- 
tonian flow defined in fl5.2l) . Moreover t„ = ti . . .t„ and /^*^^dt„ denotes the integral over he 
simplex < < < ■ ■ ■ < ti < t. Finally, E^r stands for the expectation with respect to 
the A^-particle density and 

Tj^\t) = SN{-t)Tj^^SN{t). (7.9) 
Therefore, the objects we have to investigate in the limit N oo are: 

^(Z'^^t) =$^(-1)"^ I dt„r/,(Zj;t,r,r„,t„,Z^), 

„>0 r=0 r„: rj>0 "^''^ 

\L„\=k-r 

(7.10) 

(for any configuration Zjy, typical with respect to /g'^'*), where r/j is given by: 

r/,-(Z;;t,r,r„,t„,Z^) = p2'-Ti;")(t„)4^"-)(t„_i) . . .Tp\tr)iiM,AZ-\ZN{t)). 

(7.11) 

Note that: 



,(0), 



z/r^(Z;;t) = /i^,, (z;|Z;v(t)). 



(7.12) 

us we are lead to coj 
r/i(^i;t,r,r„,t„, Zn) = P2^'4^")(t„)4^"-Ht„„i) . . .4^^)(ti)/i;v(^ll^iv(t)), 



We start by analyzing the behavior of z/j^'' in the cases j = 1,2, thus we are lead to consider: 



(7.13) 

and 

m{z[,z',;t,r,r^,t^, Z^) = V'^T^^"\QTl^"''\t^_,) . . .Tp\t^)iiM,2 {Z',\ZN{t)) . 

(7.14) 

It is useful to stress that the operators T^^\tj) (j = l,...,n) and P^*" act as suitable 
distributional derivatives with respect to the variables Z^. To evaluate ?7i, let us first analyze 
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the action of tI^\t). By (17. 9p and (13.181) . for any function G = G{Zn), we have: 
(4'-) (r) g) (Z^) = Sn (-r) (rl,^) + R%^) {S^ (r) G) (Z^) = 

]^ ^ gr+l Qr+1 

+ nY1 Y1 C^^mSn (-r) |, , - a;,) ■ , |,^| (^^ (r) G) (Z^v) • 

\kx\ + \k2\=r+l 

(7.15) 

Note that the derivatives involved here are done with respect to the variables at time t = 0. 

Denoting by any derivative of order r with respect to a variable Zj at time t = 0, we 
observe that: 

S^{-t)D:p{Z^) = {Z^{t)) = Dl^it) {S^{-t)G) {Z^), (7.16) 

where, by Dl,(t), we denote the same derivative of order r with respect to the variable Zj{t). 
Then, by ^7J6\f and (TTTSD : 



(4'-) (r) g) (Z;v) = (-r) (fi^) + i?!;)) Sj, (r) G (Z^) 



|fci| + |fc2|=r+l 

(7.17) 

Therefore, in computing the action of T^\t), we have to consider derivatives with respect to 
the variables at time r. As a consequence, we have to deal with a complicated function of the 
configuration Z^v which, however, we do not need to make explicit, as we shall see in a moment. 

On the basis of the previous considerations, we compute the time derivative of rji by 
applying the operators V'^^T^"'\tn)T^"~^\tn-i) ■ ■ ■ Tj^^\ti) to the Vlasov equation: 

{dt + v[ ■ V^'J fiNit) = (V,. * fiNit)) ■ V„i/i^(t). (7.18) 

In doing this we have to compute 

V'^T^^-\QT^^-'\t^,,) . ..T^^^\t,)fi^{z[\Z^{t))fi^{z',\Z,,{t)). (7.19) 

Now we select the contribution in which each T^p\t() and P^*" apply either on iJLN{z[\Zis[{t)) 
or to iiN{z2\Zj<4{t)). The other contribution involves terms in which are present products of 
derivatives with respect to the same variable. By Proposition 5.1 and Proposition 5.2 we expect 
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those terms to be negligible (in the C^-weak sense) in the limit iV — > oo. Therefore we obtain 
the following equation: 

{dt + v[ ■V^'^)r]i{z[,t,r, 
+ ^ 5Z {^x',(l)*m{-,t,e,rj,tj,ZN)) ■V^'^r]i{z[,t,r - e,rj^\j,tj^\j,ZN) + 

0<e<r 0<m<n ICl„: 
\I\=m, 
0<\rj\+e<k 

(7.20) 

where Ej^ is an error term which will be proven to be negligible in the limit — > cxd in Appendix 
C. In f l7.20p we used the notations: 

In = {1,2, ... , n}, I is any subset of /„, = {rj}j<zi, tj = {tj}jei- 

(7.21) 

Next, we compute the time derivative of i'[''\t). We have: 
(9, + t;i-V4)z/f) = 



E / / dt,... dt„r/i(z;;t,r,r„,t„,Z^) 

„~>n ..-n i„ I. JO Jo Jo 



n>0 r=0 |rj 



+ 

ti=t 



\rj=k~r 



+ E / dt„(9i + t;l-V,.)^i(^;;t,r,r„,t„,Z^). (7.22) 

n>0 r=0 |r,J: ^ """^ 

r/>0 
\Ln\=k-r 

In evaluating the first term on the right hand side of fl7.22p , we are lead to consider t]i evaluated 
in t = ti. Thus, according to the expression of rji (see (17.131) ). we have to deal with: 

Ti;'\t)iJM{z[\Z^{t)) = ^^(-t)4^^Viv(^il^N). (7.23) 

Therefore: 

= i-iy^^Cr, J dx'2 dv'^ V(2;'i - X'2) ■ Dll+'flN{x[,v[\ZN{t))Mx2,v'2\ZN{t)), 

(7.24) 

where the term involving off-diagonal derivatives, namely -R^^^^ (see f l3.20p ). disappears because 
both the derivatives and the empirical distribution are evaluated at time t. Hence we compute 
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7^1 in t = ti and, inserting it in the first term of tlie riglit liand side of (17.221) . we obtain: 

n>0 r-=0 r„; r,>0 "^0 "^0 "^0 

|r„|=fc-r 

= Y. (-l)^'^'cn / dv', D;;+V(a^'i " A) ■ D^j+Vf -'^^^(x;, x'^, v'^; t). 

0<ri<fc 
ri even 

(7.25) 

Let us come back now to equation (17.221) . It is useful to observe that: 

/ dt^Y. = [ / (7-26) 

J ord j^j . J ord J ord 



ICI„ 
\I\=m 



Then, putting together (I7.22p . (17.251) . (17.201) and (17.261) . we obtain the following equation for 
+ J2 (-ir^'c,, f dx',dv',D:)+'^ix[-x',)-Dlf'ut'^\x',,v[,x',,v',-t) + 

0<ri<k 
r'l even 

+ E * '^i'^ W) • ^<''t'\t) + El, (7.27) 

0<i<k 

with initial datum given by: 

u[^\x[,v[;t)\t=o = Viiiz[;0,k,ro,tQ,ZM)=V''fiNiz[\ZN). (7.28) 
Here Ejj- arises from E^^ (see (17.201) ). Now, we want to prove that: 

u[''\t) f''\t), as ^ oo, - weakly, (7.29) 

and 

i/f ^(t) ^ /f ^(t), as N ^oo, - weakly, (7.30) 

for any configuration Zj^ such that Hn ^ fo'^ in the weak sense of probability measure (namely, 
for any Z^ typical with respect to f^^)- As a consequence, reminding that h''^\t) and y'^\t) 
are equal to uj^^\{t) and uj^^\{t) respectively, a.e. with respect to W^\, (I7.29P and (I7.30p are 
equivalent to: 

^ f^\t). as iV ^ oo, - weakly, (7.31) 



and 



(k) 



(k) 



(t) ^ ft\t), as N -^oo, - weakly. (7.32) 
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As we already remarked, the C^-weak convergence implies the convergence in S', therefore, 
(I7.3ip and (17.321) imply the convergence of J^\{t) to in S'{R^ X M3) and of Jji\{t) to 



/f (t) in S'( 



7.1. One and two-particle convergence. In evaluating the behavior of I'fit) when oo, 
we note that it solves the initial value problem fl7.27l) - fl7.28l) for which we want to use Proposi- 
tion 6.2. First, however, we have to verify the assumptions. The first one, namely i), is verified 
as follows by the considerations developed in Section 5. 

Now, we have to check that assumption ii) is satisfied, namely, we have to prove that 



Vmi, U2 in C^iW X W), 

there exists a constant C = C{ui,U2) > 0, independent of A^, such that: 



Ml * {U2 i^i^\t) 



< C for any t. 



(7.33) 



We have: 

1 1 Ml * [u2 V?{t) 
k 



sup 

xi ,v'. 



dydw ui{x[ -y,v[- w)u2{y,w)iy[''\y,w;t) 



< 



< 



EE E 



t 

ord 



=k—T 



sup 



j dydw ui{x[ -y,v[ - w)u2{y,w)r]i{y,w]t]r,r^,t^, Zn) 



EE E 

\rj=k-r 



ord 



sup 

x[,v[ 



dydw {ui{x[ -y,v[- w)u2{y, w)) V'^''Tl^"\tn) . . . Tp\ti)^N{y, w\ZN{t)) 



E 

|r„l=fc-r- 



dt. 



ord 



sup 

x[,v[ 



dydw g{x[,v[,y,w)V''-Tj:^-\tn) . . .Tp\h)fiN{y,w\ZNit)) 



(7.34) 



where we used the notation g{x[,v[,y,w) := Mi(a;'i — y,v[ — w)u2{y,w) and, clearly, we have 



g{x',,v[,;-) e Q 



for any x'l and v[ and g{-, ■, y, w) G Q 



for any y and w. By 



some estimates which will be proven in Appendix C (see Lemma C.2), we are guaranteed that, 
applying the operator V^''T^^''\tn) ■ ■ ■ T^n'\^i) the empirical measure /XAr(t) and integrating 
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versus a function in C^(M^ x M^) we obtain a quantity uniformly bounded in N. This feature, 
by virtue of the good properties of the function g ensures that fl7.34p is finite. 

Let us now look at the initial datum for i'[''\t), in order to verify assumption iii). 

From (I72H1) we know that z/f^(0) = "D^Vw e S'(M.^ x R^). As regard to its limiting 
behavior, we find that: 



n=l /C/jv Sj-.jGl jel 
\I\=n l<sj<k 

where 1^ = {!,..., A^}. For our convenience, we have written the action of the operator "D^*^ 
in a equivalent and slightly different way from that we used in f l7.3p . 

We realize that the only surviving term in the sum (17.351) is that with n = 1. Hence: 

N N 

''i(t)Lo = E = ^ E ^oWi - = ^gVtv. (7.36) 
Therefore we can conclude, by using the mean-field limit: 

= (D^V/i^) - {dI'uJ!,'^) = {u,Dl'f!,'^) = (ujt^) , as AT ^ oo, 
V u in (M^ X R^) . (7.37) 

Thus, /o*^'' plays the role of 70 in Proposition 6.2 and it is in (M^ x M^) because Jq^^ G 

s{m:^ X r3). 

We conclude the convergence proof (for the one and two-particle functions) by induction. 
For = we know that, for any configuration Zn which is typical with respect to /g^'*, we 
have: 

u[''\t)=ixj,{t)^f(^\t), as A^-^00, (7.38) 

in the weak sense of probability measures, and, as a consequence, the convergence holds — 
weakly. Moreover 

4'\t) = Mt) ® Mt) /f (t) = f^'\t) ® f^'\t), as AT ^ 00, (7.39) 

in the weak sense of probability measures, and, as a consequence, the convergence holds — 
weakly. 

We make the following inductive assumptions for all h < k: 

u[''\t) j<''\t), as N -^00, - weakly, (7.40) 
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for any configuration Z^- which is typical with respect to flf \ and 

^W(t) ^ /f = ^ ^(.)(^)^(/.-.)(^)^ as iV ^ oo, C - weakly, (7.41) 

0<g</i 

for any configuration which is typical with respect to /o^"*. 

Now we want to prove that (17.401) and (I7.4ip hold also for h = k. 
Thanks to (I7.40p . we can affirm that: 

o<e<k o<e<k o<e<k 

- weakly, (7.42) 
and, thanks to (I7.4ip . have: 

0<ri<fc 
ri even 

i - weakly 

— <k 

ren 



0<ri<fc 
ri even 



0<ri<A; 0<i3<fc-ri 
ri even 



j2 j2 T^'''^f^'-'''-'\t). 



0<ri<k 0<q<k-ri 
ri even 



(7.43) 

At the end, putting together (I7.42p and (I7.43p . we find that the sum of the source terms 
in equation fl7.27p converges C^-weakly to: 

0<e<k 0<ri<fc 
0<5<fc— ri 

which plays the role of G in Proposition 6.2 and it is easy to check that it is in C° (L^(R'^ x M^), R+). 
Therefore, we can apply Proposition 6.2 claiming that, for any typical configuration Zn with 
respect to /o°'', i'[''\t) converges C^-weakly to the solution of the problem (I6.10p . Looking at 
(I2.15P and (I2.17p . we realize that we obtained the equation satisfied by f^''\t). 

In order to "close" the recurrence procedure, it remains to show the two-particle conver- 
gence at order k. It follows from the one-particle analysis and from the following computation 
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(see (Ol): 

r]2{z[,z'2;t,r,r^,t^,ZN) = 

= ^ Y Vi{z[]tJ,rj,tj,ZN)vi{z2,t,k - e,rj^\j,tj^\j,ZN) + R%, (7.45) 

0<e<k 0<m<n I:ICI„ 
\I\=m 

where i?^ is a remainder arising from the action of the operator P^'^T^"-'(t„) . . . Tj^^\ti) on a 
product of two empirical measures fiN{t)- In Appendix C we will see that it is vanishing in the 
limit. As a consequence, 1^2''^ (see (17.101) for j = 2) is such that: 

ui'\t)= J2 ^['\t)ut'\t) + o{l), (7.46) 

0<q<k 

in the limit N ^ 00. Therefore, from the inductive assumption f l7.40p and from the one-particle 
convergence at order k, we conclude that: 

4'\t) ^ E f^'\t)f'^'''\t) = (^)' as N^oo, Cr - weakly, (7.47) 

Q<q<k 

for any configuration Zn which is typical with respect to Thus, we have just proven the 
convergence of u^^^- in the cases j = 1, j = 2. 

7.2. j-particle convergence. As for j = 2, the j-particle convergence can be reduced by 
the one-particle control. Indeed by f l7.10p and (17.111) we have: 

si-.-Sj m=l 

0<Sm^k 
Em ^rn = k 

with i?-^ when 00. 

Again the error term i?]^ arises from the presence of products of derivatives with respect to 
the same variable. In conclusion, the result we proved for i'[''\t), together with the estimates 
proven in Appendix C, is sufficient to guarantee the C^-weak convergence of i'j''\t) to ff^\t) 
for any j (for any typical configuration Z^ with respect to /o°^), and, as a consequence, the 
C^-weak convergence of uj^^^-{t) is fj''\t), for any j. 

The final step is to realize that this convergence does imply that for the coefficients W^j{t), 
namely what is established by Theorem 6.1. 



(7.48) 
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First of all, we observe that, for any test function u we have: 

dZNWi!;\z^;t)uizi) = 

K3JVxR3]V 

1 ^ 

dZ^Wj^\z^;t)-y^u{zi] 

dZNWji'\ZN;t) {u,fiN) = [u,uj']^{{t) 

(7.49) 

where we made use of the symmetry of the coefficient wl^\Zj^; t) with respect to any permuta- 
tion of the variables (the computation is the same we did in Section 5 for wl}\{t)). From (\TA9\\ . 
we can see that W^\{t) and uj^^\{t) are equal as distributions in S' {M.^ x R'^) (in particular, 
we can choose test functions belonging to (M'^ x M'^)), then the convergence of W^\{t) is 
proven. Moreover, for j > 2, a straightforward computation shows that, by fixing an index j, 
we have 

m,^,\_N{N-l)...(N-3 + l)f_^ wWm\ , ^.<J 



-J. = ' "'^ - ' K^m) + (7.50) 

where Cj^-j < oo provided that (^Uj,W^^j(t)^ is uniformly bounded for each j < j. Then, to 
conclude the proof of Theorem 6.1, it is enough to use a recurrence argument. 
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Appendix A 



Proof of Proposition 5.1 



To avoid inessential notational complications, we deal with the one-dimensional case. 
By the Newton equations, we have: 



dVr 
dVr 



5„.t+ / ds{t- s)— Y^d.,F 

Jo -'^ ■ / • 



dxi{s) dxj{s] 
dvr 



dvr 



Xj{s)) 



dxi{s) dxj{s] 



dvr 



dvr 



where: 



-V,(; 



(A.l) 
(A.2) 

(A.3) 



is the force associated with the potential (f). 

Let us analyze in detail the derivative of Xi{t). From (lA.ip . we get: 



max 

i,r 
t<T 



dxi{t) 



dVr 



< c. 



(A.4) 



Inserting this estimate again in (lA.ll) . we realize that we can obtain a better bound for 
in the case r i (see namely: 



dxi{t) 



dVr 



< C / dsit - s) 



dxi(s) 



+C / ds(t-s) 



dVr 

dxAs 



+ 



(A.5) 



dVr 



-C / ds(t-.)-5^9.F 

J ■ / ■ 



+ 



[Xi{s) - Xj{s)) 



dxj{s) 



dVr 



Hence, by virtue of the Gronwall lemma, we find: 

dxiit) 



max 

t<T 



dVr 



< 



C_ 
N' 



(A.6) 



(A.7) 



By (IA.2p . we find that the same estimate holds for the derivative of Vi{t) with respect to Vr- 
Analogous estimates hold for the derivatives with respect to the initial positions (see also |17j). 
Therefore the claim of Proposition 5.1 is proven for derivatives of order one. 
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Now, let us consider a sequence / := (ji, . . . ,jk) of possibly repeated indices. We show 



that: 



N 



N ^ 



i=l 



d^Xi(t) 



dvj, ■ ■ ■ dvj^ 



< 



C 

Ndk ' 



(A. 



where dk is the number of different indices in the sequence ji, . . . We know that flA.SP is 
verified for = 1 (it follows directly by (IA.4p and flA.7P ). thus we prove (lA.SP by induction on 
k. Denoting by: 

m ■■= . , (A.9) 



dvj, . . . dvj^ 



estimate (lA.Sp can be rewritten as: 



1 ^ 

-5^|D(J)x.(t)| 

1=1 



< 



C 



(A.IO) 



By (lA.ip we derive the following estimate for D{I)xi{t): 

\D{I)xM< d.(t-s)-X;i^a)(^.(^)-a:,(s))|+M,(t), 



(A.ll) 



where the term Mi{t) can be computed from (lA.ll) according to the Leibniz rule. Let Vn 
{/i, . . . , /„} be a partition of the set / of cardinality n, with 2 < n < k, then we have: 

N k 

j^i n=2 
k - ^ , 

' n [D{H) {x,{s) - x,{s))] 



t N k 

M,{t)< / ds(t-s)-5^5^5^C(P, 

"^0 j^i n=2 Vr. 

pt k N 

"^0 n=2 P„ j=l 



< 



(A.12) 



where D{H) := Y\.heH C{Vn) are coefficients depending on the partition Vn and on 
suitable derivatives of F. By ( lA.llI) . it follows that: 

-J2\DilMt)\< / ds{t~s)-J2\DilMs)\ + Mit), (A.13) 
i=i i=i 



where M{t) = ^ J2i=i ^i(^) and, by ([02]) . we have: 

( fe N N 

"^0 n=2 V„ i=l J=l 



(A. 14) 
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We observe that: 



1 ^ 
1 ^ 



HeVn 



1 ^ 



i=l HdVn 



QcPn V i=i QeS / y i=i JeVn\Q 



where Q is any subpartition of Vn and C(Q) are coefficients depending on Q. 
We assume that the estimate (lA.lOp holds for any m < k — 1, namely: 

1 ^ C 

- J2 \D{M)xi{t)\ < j^, for any McI s.t. \M\ = m < k - 1, 



(A.15) 



(A.16) 



i=l 



where dm is the number of different indices in the sequence M. 

Indeed, if we consider a partition Vn of cardinality n > 2, we are guaranteed that \M\ < k — 1 
for each M &Vn- Then, by noting that: 

1 ^ 1 ^ 

]^ E n \^(H)x^it)\ < n ]^ E \D(H)xiit)\ , V subpartition H C p„ (A.17) 



1=1 Hew 



N 



we can apply the inductive hypotheses flA.16P to estimate the derivatives of Xi{s) and Xj{s) 
appearing in flA.151) . Thus, we obtain: 



n WH) 



Xi{s)\ < 



1 ^ 

n i^Ei^(^)-^' 



i=l H£Vn 



N 

HeVn i = l 



s)\ < 



< 



n 



C 



c c 

< 



(A. 18) 



where dh is the number of different indices in the sequence H and we used that J^HePn dk- 
In a similar way, we find 



1 ^ 

^Eni^(^)^^(^)i^ n 



c 



(A.19) 
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where dq is the number of different indices in the sequence Q. 
Moreover, we have: 



1 ^ 

n i^(^)-.(^)i< n 



j=i Her„ 



C 



C C 
< 



(A.20) 



and 



C 



(A.21) 



n i^(^)-.(^)i< n 

i=i JeP„\Q JeVn\Q 

where dj is the number of different indices in the sequence J. Then, putting together (1A.19I) 
and dA^Il), we find: 

E ^(2) f^E n \D{QMs)\] ( ^E n ) < 



Ec^(s)n n 



C 



< 



]S[dq+dj 

QeQJeVn\Q 

s E^ran n 

QcPn QeQJeV„\Q 



(A.22) 



In the end, we have just proven that each term in (lA.lSp is bounded by Therefore, by 

using this estimate in (1A.14I) . we find: 



M(t) < 



C 



(A.23) 



By ([03]) and f[03|l . it follows that: 

-Ei^aKWi< I d.(t-s)^Ei^(^)^^(^) 



c 



Therefore, by using the Gronwall lemma, we find: 



1 ^ 

-Ei^w^.wi 

4 = 1 



< 



c 



(A.24) 



(A.25) 



. . . , Vj^, an 



As regard to the derivatives of Vi{t) with respect to some initial velocities v 
analogous estimate holds and the proof works in the same way. Furthermore, this strategy 
leads to the same estimate for the derivatives of the function J2iLi ^ii'^) with respect to some 
initial positions , . . . , Xj^. . 
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Now, thanks to the estimate we have just proven for the derivatives of the function 
F Si^i -^*(^)' ^t)le to prove the claim of Proposition 5.1. In fact, we have: 



N N N 

i=l i=l 1=1 



i=l 



c 



(A.26) 



where Del contains the different indices appearing in the sequence /. Thus, according to our 
previous notation, \D\ = dk and we denote the elements of D by ji, . . . , jd^.. Then by flA.261) 
we find: 



1 ^ 1 1 

- E \D{I)zm = ^ \D{I)z~,M + ... + - \DiI)z.^^it) 

i=l 



+ 



1 ^ 



i=l 



C 



(A.27) 



which implies 



\D{I)zm<C 



+ 



or 



c 



\D{I)z,{t)\< 



(A.28) 



(A.29) 



where d!)^ is the number of different indices in the sequence I which are also different from i. 



□ 



Appendix B 



Proof of Proposition 6.1 

Let f/h(t, s) be the two parameters semigroup solution of the linear problem: 

{dt + v- V,) Uhit, s)7o = (V0 * h) * VMh{t, s)jo, 
Uhis,s)-fo = 7o. 



(B.l) 



The solution of ( IB.ip is obtained by carrying the initial datum 70 along the characteristic flow 

X = V, 

V = — V0 * h. 



(B.2) 
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Next, we consider the problem 

7|t=o = 7o- 

which can be reformulated in integral form: 

7(t) = Uh{t, 0)70 + [ ds Unit, s) [(V0 * ^{s)) ■ V./i(s)] . (B.4) 
The above formula can be iterated to yield the formal solution 

^(x,v;t) = Uh{t, 0)'yo{x, v) + ^ / ^'^i / ■ ■ ■ dt„ / dxi / dt;i . . . / dx„ / dvn 

Jo Jo Jo J J J J 

Uhit, ti) [V^h{x, v; ti) ■ V^0(x - xi)] 

Uh{ti,t2) [Vv^h{xi, vi; t2) ■ Vx^<P{xi - X2)] 

Uh{tn-1, tn) ["^Vn-iH^n-l, Vn-l, tn) " x„-i4>{Xn-l " X^)] 

)• (B.5) 

We remark that Uh{tk, tk+i) acts on the variables Xk, Vk with the convention that (xq, fo) = (a;, f ) 
and, furthermore, Uh is multiplicative and preserves the Lp(M^ x M^) norms (p = 1,2,..., 00). 
Under the assumptions of Proposition 6.1, the above series is bounded in L^(M^ x M?) by: 

/ . ^ (^SUp^g[o^t] ||Vi,/l(r)||^l(]g3^]g3) j ||V:;;0|lLoo(lf53) ||7o|Il1(M3xK3) , 
n>0 



(B.6) 



which is converging for each t. Now, we denote by S/i(t, s) : L-'^(]R^ x M?) L^{E? x R-'^), the 
two parameters semigroup given by the series (IB.Sp . Then, the solution 7 to the problem (16.101) 
is given by: 

7(t) = 0)70 + /" ds llhit, s)e(s), (B.7) 







and, thanks to the assumption we made on B and to the fact that the above series fIB.Sp is 
converging for any t, we are guaranteed that 7 G C° (L-'^(]R^ x R^),]R"'"). 

The regularity of 7(t) = S/i(t,0)7o follows by fIB.Sp and the fact that Uhit.ti) propa- 
gates the regularity. 

□ 



Proof of Proposition 6.2 



(B.8) 
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The proof consists of two steps. 
Step 1): 

Let be as in Proposition 6.2. Then, we show that solves the problem: 

{dt+V- V^) 7Ar = L{h)-iN + 0^, 
7Ar|t=0 = iNfl-, 

with 

e^ = e,v + i?iv, (B.9) 

and i?Ar is such that: 

Rn 0, - weakly. (B.IO) 
In proving (IB.lOp . the assumption ii) on 7^? is crucial. 

Step 2): 

By virtue of Step 1), the hypotheses we made on Vyh and Proposition 6.1, we find that: 



7^(t) = 0)7Ar,o + I ds s)Q'j,{s). (B.ll) 

^0 



Then, reminding that: 

o h{t) e C^{M? X R3) for any t, 

o the flow propagates the regularity, 

o i?^ ^ 0, - weakly, 

and by virtue of the assumptions on 7Ar o and 6^^, we can easily show that: 

7^ ^ 7, as N ^ 00, - weakly, (B.12) 

where 

7(t) = S;,(t,0)7o+ / dsT.h{t,s)Q{s). (B.13) 

Therefore, we recognize that 7 solves the problem (16.101) and, by virtue of Proposition 6.1, it 
is uniquely determined by flB.lSp and hence it is in {L^{E? x ]R^),M+). 

Proof of Step 1 ): 
We have: 

{dt + v- V^) 7Ar = L{h)'yN + ©tv + - /i)77v 

/ N (B.14) 

lN{x,v;t)\^^Q = 7jv,o(a;,t;), 

where 

Rn = Rn{x, v; t) := L{hiy — h)'jN- (B.15) 
We want to show that Rn 0, C^-weakly. According to the definition of the operator L, we 
have: 

Rn = (V,0 * {hN - h)) V,7iv + (V,.0 * 7^) V,{hN - h), (B.16) 
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thus, we have to show that 

{u, * (hN - h)) V^7jv) 0, as oo, V n G x M^), (B.17) 

and 

{u, (V^(/) * 77v) V^(/i7v - /i)) ^ 0, as ^ oo, V u G C^{R^ x R^). (B.18) 
We show only (IB.lSp in detail because (1B.17P will follow the same line. We have: 



{u, (V^0 * 7Ar) V^{hN ~ ^)) = y ^^^'^ J ^y^'^ ^) - y)lN{y, w; t) ■ 

Vi,(/lAr(x, V] t) - h{x, v; t)) = 

{h^ix, v; t) — h{x, v; t)) = 
dxdv J dydwVvu{x,v) (Vx0 * 7Ar) {x,v;t){h — hN)ix,v;t). 

(B.19) 

Setting 

Cn{x,v) ■.= V^u{x,v) j dydwV^(j){x - y)'jNiy,w;t), (B.20) 
we can write (IB. 191) as: 

{u, (V^0* 7Ar) Vy{hN -h)) = j dxdv CNix,v){h{x,v;t) - hN{x,v;t)) = 
dxdv / dx'dv' {(Nix,v) — C,n{x' ,v')) Pn{x,v; x' ,v';t), 



(B.21) 

where Pn is a coupling of h and hpf, namely a probability density in x with marginals 
given by h and h^- Now we observe that: 

V^,^CAf(x, v) := J dydwVx,v [Vj,m(x, t;)V^0(a; - y)] 7Ar(y, t), 

(B.22) 

and, thanks to the assumption ii) we made on 77V, we know that there exists a constant 
C = C{u, (p) > such that: 



sup |V^,^,CAr(a;,t;)| = || VCjv ||ioo(B;3x]R3-, <C < +00. 



x,v 



(B.23) 
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Therefore, coming back to ( 1B.21I) . we find: 



|(M,(V,0*7iv)V,(/iiv-/i))| ^ y y dz'\Cj,{z)-C,,{z')\PN{z;z';t) 

< j j dz'C\z- z'\Pn{z; z';t). 



(B.24) 



wfiere we used the standard notation z = {x,v) and z' = {x',v'). Then, taking in ( ]B.24p the 
infimum over all couplings between h and h^, we obtain that: 

\{u, (V,0 * 77v) V,(/i^ -h))\< CWih^, h), (B.25) 

where, as in Section 5, W denotes the Wasserstein distance. But we know that the right hand 
side of (1B.25I) goes to zero because of the assumption i), then we have just proven that: 



\{u, (V,0 * 7^) ^v{hN -h))\-^0, V M G C(R=' X 
Analogously, we can prove that 



\iu, (V,0 * {hN - h)) V,7^)| ^0, V M G C(M' X 



(B.26) 



(B.27) 



Therefore we have just proven that Rn goes to zero in the C^-weak sense and the proof of Step 
1) is done. 

Proof of Step 2): 

Thanks to Step 1) and to the assumption on V„/i, we know that 7Af(t) can be written as 
in (IB. lip . Then, for any function u in C^(]R^ x M^), we have that: 



(u, 7^(t)) = (u, S^(t, 0)7^,0) + / ds (m, S^(t, s)e'^(s)) , (B.28) 

namely 

K7Jv(t)) = ((S^(t,0))*M,7jv,o)+ / ds ((S;,(t,s))*M,0^(s)), (B.29) 







where SJj is the adjoint of S/^. We remind that the two-parameters semigroup s) prop- 
agates the regularity, provided that V^h G C*^ (R^ x M^). In particular, if acts on a 
function u which is in C^(]R'^ x M'^) and the function h{t) is supposed to be in C^(]R^ x M^) 
for any t, as it is in the assumptions of Proposition 6.2, we are clearly guaranteed that Vyh{t) 
is in C^{R^ x M^) for any t, and then, u{t) := 0)m(x, v) is also in C^(R^ x M^) for any t. 
Obviously, the same holds for Sjl^. Thus, the functions 0))*m and {T,h(t, s))* u appearing 

in flB.29p are in C^(R^ x M.^) for any t. Therefore, thanks to the assumptions we made on 77v,o 
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and Gat, and of what we know about Rn, we find that: 

((E;,(t,0))*M,77v,o)+ / ds {{Eh{t,s)Tu,e'^{s)) 

Jo 

I N -^oo 
t 



((E^(t,0))*M,7o)+ / ds ((S^(t,s))*M,e(s)) 



u,S;,(t,0)7o)+ / ds {u,T.h{t,s)Q{s)). (B.30) 







Finally, by Proposition 6.1, we know that the expression ( IB. 301) identifies properly the unique 
solution of the problem fl6.10p in [L^ {M? x M?) ,M+) and Proposition 6.2 is proven. 

□ 



Appendix C 

Lemma C.l: For each time r > 0, let us define the operator T^\t) as follows: 

Then, for each m > and for each u G C^(M'^ x M?), there exists a constant C > 0, not 
depending on N , such that: 

\)\(u,f^-j-\t^)...fp\t^)^iM{t) 

Moreover, we have: 



< C. 



(C.l) 



ii) {u,T^^-\t^)...TP>{t,)fiN{t)) < {u,T'^-\tm)...T'^"{ti)fiN{t) 



(C,2) 



Proof: 

We observe that: 



u,fj^-\t^) . ..f^^^\t^)Mt)) = fl;-\t^)f^^-^\t^_,) . ..fp\t,)U{Z^{t)), 



where: 



1 ^ 

U{ZN{t)) := iu,fiN{t)) = 



(C.3) 



(C.4) 
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We assume m > being the case m = obvious. 
By using the notations: 

5(r^,tJ:=4""^(U---4"Hti) (C.5) 

and 

S{r^,tJ:=fl^-\t^)...f^^^\t,), (C.6) 
we have (see the first term in the right hand side of (17.151) ): 

jl---jm ll ■••Ijn 



Dl^^'<P{x,,{t,)-xiAti)) ■ D:]^\t,)U{Z^{t)), (C.7) 



C depending on r^. By setting: 



1 ^ 

<l>,„(Z^(t„)) := ^ 5^ D:-^'<P{xM - ^iM) (C.8) 



N 

ln = l 



V n = 1,2, ... ,Tn 



(IC.7P can be rewritten as 



<l>,_,(Z^(t„_i))-Z};;"-+i(t„_0 



^,,{Z^{t,))-Dl^^\t,)U{Z^{t)). (C.9) 

We observe that, thanks to the smoothness of the potential (f), $j„ (for each n) is a uniformly 
bounded function of the configuration Z^, together with its derivatives. 

Performing the derivatives in (lC.9p . we realize that S{v_^,t^U{Zi^{t)) is a linear combination 
of terms of the following type: 



D:-HU...D5^(t2)$,,(^^(ti)) 

D:2"'{tJ<l>,^_,{Z^{t^.,)), (C.IO) 
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with the constraint 



0-1,1 = ri + 1 

0-2,1 + 02,2 = r2 + 1 



(C.ll) 



''171,1 



+ Om,2 + ■ ■ ■ + Ctri 



+ 1. 



For a fixed sequence we have to compensate the divergence arising from the sum 
^^.^ J , which is 0{N"^), by the decay of the derivatives as given by Proposition 5.1 and 
Proposition 5.2. Indeed we have: 

D:2'{tJ...D:^^^{t,)D:^^^{t,)U{Z^{t))\ < ^, (C.12) 

where d is the number of different indices in the sequence ji,j2, ■ ■ ■ ,jm for which a^,!, . . . , 02,1, ^i^i 
are strictly positive. Note that the fact that the derivatives are not computed at time t = 



but at different times ti, t2 



,tm, does not change the estimate in an essential way. 



An analogous estimate holds when we replace U by some namely 



< 



C 



which are also different 



where d^-i is the number of different indices in the sequence j^, . . . , 
from and from which am,k, ■ ■ ■ , Ok,k are strictly positive. 

As regard to the term in the sum • . in which all the indices are different (which is 
the only one of size 0{N"^)), the constraints (IC.llI) together with estimates (]C.12p and (IC.ISP 
ensure that the product of derivatives on the right hand side of (IC.lOP is bounded by 



Thus this term is of order one. Now for each s 



,m — l consider the 



terms in the 



s!(m,— s)! 

sum Ylji j„ which s indices are equal. The sum is bounded by N"^~^. On the other hand, 
the constraints fIC.lip together with flC.12p and flC.lSp ensure that the product of derivatives 
on the right hand side of fIC.lOp is bounded by 1/N'^~'^ . Thus even these terms are of size one 
and i) is proven. 

To prove ii) we observe that: 



S{v„,tjU{Z^{t)) - S{v^AjU{Z^{t)) (C.14) 
can be expanded as in ( ]C.7p and ( IClOj) . However now we have an extra derivative, arising 



from the definition of (see (I3.20p ). which yields an additional We omit the details of 
the proof which follows the same line of i) . □ 



In the same way we can also prove the following 

Lemma C.2: For each m > 0, k > and u G C^(M^ x M^), there exists a constant C > 0, not 
depending on N, such that: 

\V''S{r^,tjU{ZMm<C. (C.15) 
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where U{ZN{t)) is defined as in ( [(7.^[ j. 



Proof: 

First we look at the case m > 0. Reminding the structure of the operator V"^^ (see (17. 3p ). 
we are led to consider the term D^*^S'(r^, t^)f/(Z7v(t)). We remind that D^q^ is a derivation 
operator with respect to the variable Zj that acts as specified by fl6.8p . By the expansion (IC.lOp 
we readily arrive to the bound: 



£'GjS(r„,t„)f/(Zjv(()) 



- N 



(C.16) 



Indeed by applying Dj^- to (IC.lOp either j ^ (ji ■ ■ . jm) so that we gain by the extra 
derivative, or j G (ji . . . jm) so that we reduce the sum 'Y^j^ by a factor More generally, 
by the same argument we find: 



< 



c 



(C.17) 



where n = 

Finally by writing the action of the operator P^'^ as in fl7.35p . we obtain 



N 



N 



n=l ^ ' Sl---Sn 

l<Sj<fc 



^ n\(N - n)\ iV" 

n=l ^ ' 



c 



< 1 + ^ 1 <c. 



N 



(C.18) 



B, C being positive constants not depending on N. Again V^''S{r^, t^)U{ZN(t)) is the leading 
term of X'^'^5'(r^, t^)f/(Zjv(t)) for the same reasons we discussed in Lemma C.l. 

If m = 0, the estimates flU.16p and flC.17p follow directly by Proposition 5.2. Thus, even 
in this case, the proof is concluded by ( IC.lSp . 

□ 



The fact that the error term E]^ (see (\T20^ ) and hence E% (see (17:271) ) are C^^-weakly 
vanishing when — »• oo is an immediate consequence of the following 



Lemma C.3: Let r j and tj be defined as in Section 7, for any J (Z In with J„ = {1, 2, . . . , n}. 
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For any r > we have: 

V''-S{r^,tJfiN{z[\ZN{t))fiN{z',\ZN{t)) = 



r,N 



0<£<rO<rn<n /C/„ 
\I\=m 



where 



er,N — as N ^ oo — weakly. 



(C.19) 
(C.20) 



Proof: 

It is enough to prove (1C.19I) and (1C.20I) replacing each streak S with the corresponding S, 
being the difference S — S neghgible in the hmit. 

We start by assuming r = 0. In that case, testing the left hand side of (jC.191) against a 
product of two test functions Ui,U2, we are led to consider: 



S{r„,tjU,{Z^{t))U2{Z^{t)) 



for which we can apply the expansion flC.7p . 

Proceeding as in the proof of Lemma C.l (see fIC.lOp ). we have to consider: 

D^-n^™) • • • D^,l^it2)D'^^^^ih)U,iZ^it))U2iZ^it)), 

where ai^i = ri + 1 > 0. Now any contribution of the form 

it,)U,iZ^it))Dl (ti)f/2(^^(t)), 



(C.21) 



(C.22) 



(C.23) 



with a>0, P>0, a + P = ai i is O {^), therefore it is negligible in the limit. The same 
argument applies to Dvf^{tk) whenever a^^i > . This means that each derivative appearing in 

5* either applies to HN{,z[\ZN{t)) or to /iAr(z2l^Af(^)) up to an error eo,Ar vanishing in the limit. 
This is exactly what (laToD and (1020]) say for r = 0. 

For r > we have to apply V'^'^ to (1C.19I) (replacing S by S) with r = 0. Clearly V'^^eo^N 
vanishes in the limit. Moreover: 



D 



2s i 



G,j 



5(r„t,)f/i(Z^(t))5(r,„\„t,„\,)t/2(Z;v(t)) 
D'^^^S{r,,tj)U,iZ^it))) Sirj^\„tj^\j)U2iZ^it)) + 

By simple algebraic manipulation we finally arrive to (ICIOP and (lC.20p . 



+S{rJ,tJ)U^{Z^{t)) {D'j;^S{rj^^„t,^^j)U,{Z^it)))+0{ ^ 



(C.24) 
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